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It is shown that a QCD based nuclear absorption model, with few parameters fixed to reproduce 
experimental J/ip yield in 200 GeV pp/pA and 450 GeV pA collisions can explain the preliminary 
PHENIX data on the centrality dependence of J/ip suppression in Cu+Cu collisions at RHIC energy, 
Y / sjvjv=200 GeV. However, the model does not give satisfactory description to the preliminary 
PHENIX data on the centrality dependence of J/ip suppression in Au+Au collisions. The analysis 
suggest that in Au+Au collisions, J/ip are suppressed in a medium unlike the medium produced in 
SPS energy nuclear collisions or in RHIC energy Cu+Cu collisions. 
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PACS numbers: PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 



Lattice QCD predicts that under certain conditions 
(sufficiently high energy density and temperature), or- 
dinary hadronic matter (where quarks and gluons are 
confined), can undergo a phase transition to a decon- 
fined matter, commonly known as Quark Gluon Plasma 
(QGP). Nuclear physicists are trying to produce and de- 
tect this new phase of matter at RHIC, BNL. J/ip sup- 
pression is recognized as one of the promising signal of 
the deconfinement phase transition. Due to screening 
of color force, binding of a cc pair into a J/ip meson 
will be hindered, leading to the so called J/ip suppres- 
sion in heavy ion collisions However, J/ip's are ab- 
sorbed in nuclear collisions also and prior to NA50 158 
AGeV Pb+Pb collisions |2|, all the experimental data 
on J/tp suppression, are explained with nuclear absorp- 
tion only. NA50 collaboration measured centrality de- 
pendence of J/ip suppression in 158 AGeV Pb+Pb colli- 
sions. Data gave the first indication of 'anomalous' mech- 
anism of charmonium suppression which goes beyond the 
conventional nuclear absorption. The data generated lot 
of excitement as it was believed to give first indication 
of QGP formation. Later, the data were explained in 
a variety of models .with or without the assumption of 
QGP g 1 i, i 0, i, ©, M EL G1 El- More recently, 
NA60 collaboration measured the centrality dependence 
of charmonium suppression in 158 AGeV In+In collisions 
OGHIlil]- I n In+In collisions also, one observes anoma- 
lous suppression, which is beyond the conventional nu- 
clear absorption. 

In recent Au+Au collisions at RHIC, one observe 
a dramatic suppression of hadrons with high momen- 
tum,transverse to beam direction (high pr suppression) 
0, E, EH H3|- This has been interpreted as evidence 
for the creation of high density, color opaque medium of 
deconfined quarks and gluons [2l[. It is expected that 
high density, color opaque medium will leave its imprint 
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on J/ip production. At RHIC energy, it has been ar- 
gued that rather than suppression, charmoniums will be 
enhanced [12, IH| . Due to large initial energy, large num- 
ber of cc pairs will be produced in initial hard scatter- 
ings. Recombination of cc can occur enhancing the char- 
monium production. PHENIX collaboration have mea- 
sured the centrality dependence of J/ip invariant yield 
in Au+Au collisions at RHIC energy, ^/snn— 200 GeV 
[H, d^]. More recently, with improved statistics, they 
have measured J/ip's in Cu+Cu and in Au+Au collisions. 
Preliminary results for the centrality dependence of nu- 
clear modification factor (Raa) and mean square trans- 
verse momentum for J/ip suppression in Cu+Cu and in 
Au+Au collisions are available [IE US- PHENIX data 
on J/ip production in Au+Au/Cu+Cu collisions, are not 
consistent with models which predict J/ip enhancement 
[22l l23l ]. It was also seen that various models, e.g. co- 
mover model H, statistical coalescence model or the 
kinetic model [5| fail to explain the PHENIX (prelimi- 
nary) data on the nuclear modification factor for J/ip in 
Cu+Cu and in Au+Au collisions. The data are also not 
explained in normal nuclear absorption model [28| . 

We have developed a QCD based nuclear absorption 
model to explain the anomalous J/ip suppression in 158 
AGeV Pb+Pb collisions [1, [i~3l ]. Unlike in conventional 
nuclear absorption model, in the QCD based nuclear ab- 
sorption model, the cc pair interact with the medium 
and gain relative 4-square momentum. Some of the pairs 
can gain enough 4-square momentum to cross the thresh- 
old for open charm meson, reducing the J/ip yield. The 
parameters of the model were fixed to reproduce J/ip 
yield in pp and pA collisions. The model give consistent 
description of the centrality dependence of the J/ip sup- 
pression and px broadening in 158 AGeV Pb+Pb col- 
lisions and in 200 AGeV S+U collisions [ll|- In the 
present paper we have tested the model against the pre- 
liminary PHENIX data on J/ip suppression in Cu+Cu 
and Au+Au collisions at RHIC energy, ^/snn—^OO GeV. 
Centrality dependence of J/ip suppression, in Cu+Cu 
collisions, is well explained in the model, but the model 
fails to explain the suppression in Au+Au collisions. The 
analysis suggests that in Au+Au collisions at RHIC, 
J/ip's are suppressed in a medium unlike the medium 
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produced in S+U/Pb+Pb collisions at SPS energy or in 
Cu+Cu collisions at RHIC energy. We also apply the 
model to explain the preliminary PHENIX data on cen- 
trality dependence of px broadening for J/tp's. Within 
errors, pr broadening at RHIC seems to be consistent 
with that at SPS energy. 

The paper is organised as follows: in section II, 
we briefly describe the QCD based nuclear absorption 
model. In section III, PHNIX data on the centrality de- 
pendence of J/i/j suppression in Cu+Cu and in Au+Au 
collisions are analysed. Centrality dependence of pt 
broadening of J/i/j 's are analysed in section III. Summary 
and conclusions are drawn in section IV. 



II. QCD BASED NUCLEAR ABSORPTION 
MODEL 

In the QCD based nuclear absorption model [§, 
11 , [13, E3| , J I ip production is assumed to be a two step 
process, (a) formation of a cc pair, which is accurately 
calculable in QCD and (b) formation of a J/ip meson 
from the cc pair, a non-perturbative process, which is 
conveniently parameterized. The J/tp cross section in pp 
collisions, at center of mass energy ^fs is written as, 
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FIG. 1: Energy dependence of total J/tp cross-section in pp 
collisions. The solid line is the fit obtained to the data with 
Eq.l. 
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where b runs over all parton flavors, and Q 2 — q 2 + 
Am 2 . The K factor takes into account the higher order 
corrections. We have used the CTEQ5L parton distribu- 
tion function for <p(x, Q 2 ) (29| . The incoming parton mo- 
mentum fractions are fixed by kinematics and are; x a — 
(y^| + AQ 2 /s+x F )/2 and x b = {y/x*. + AQ 2 /s~x F )/2. 
dab^cc are the sub process cross sections and are given 
in 30]. F c5 _ > j^(q 2 ) is the transition probability that a 
cc pair with relative momentum square q 2 evolve into a 
physical J/ip meson. It is parameterized as, 
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All the energy dependence of J /if) production is con- 
tained in EqflJ As shown in Fig[TJ with KNj/^ as a 
overall normalisation, over a wide range of energy, in- 
cluding the RHIC energy, the model correctly reproduces 
the experimental J/ip cross sections in pp collisions. 

Our main interest here is to test the model against the 
preliminary PHENIX data on the centrality dependence 
of J/tp suppression in Cu+Cu and Au+Au collisions at 
RHIC energy [2^,|27|. In AA collisions, at impact param- 
eter b, number of J ftp mesons produced is calculated as, 



N%+Q>) = a J Jt J d 2 S T A (s)T B (b - s)5(L(b,s)), (3) 
where Ta,b are the nuclear thickness function, 

T A {h) = J dz PA (b,z), (4) 
For the density we use the Woods-Saxon form, 
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l + exp((r-R)/a)' 



d 3 rp(r) = A (5) 



with i?=6.38 (4.45) fm and a=0. 535(0.54) fm for Au(Cu) 
nucleus (3lj |. 

In Eql2l S(L) is the suppression factor due to pas- 
sage of J ftp through a length L in nuclear environment. 
As mentioned earlier, in the QCD based nuclear absorp- 
tion model, J/ip J s are suppressed due to gain in rela- 
tive 4-square momentum of a cc pair. In a nucleon- 
nucleus/nucleus-nucleus collision, the produced cc pairs 
interact with the nuclear medium before they exit. Inter- 
action of a cc pair with the nuclear environment increases 
the square of the relative momentum between the pair. 
As a result, some of the cc pairs can gain enough relative 
square momentum to cross the threshold to become an 
open charm meson. Consequently, the cross section for 
J/tp production is reduced in comparison with nucleon- 
nucleon cross section. If the J/tp meson travel a distance 
L, q 2 in the transition probability is replaced to, 
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z 2 L, 



(6) 



e: 2 being the relative square momentum gain per unit 
length. The length L(h, s) that the J/tp meson will tra- 
verse is obtained as, 



L(b,s) =n(b,s)/2 Po (7) 
where n(b,s) is the transverse density, 

n(b, s) = T A (s)[l - e - CT ««T B (b-s)] + ^ B ] ( 8 ) 

J/tp suppression in the model is governed by the pa- 
rameter e 2 and L (see EqJSJ) . The length L is a geometric 
term. It has weak energy dependence from the energy 
dependence of a n n, the inelastic NN cross-section. En- 
ergy dependence of J/tp suppression will reflect, most on 
the parameter e 2 . NA50 data on J/tp production in 450 
GeV pp/pA collisions and in 200 GeV pA collisions are 
well fitted with a common square momentum gain factor, 
e 2 =0.187 GeV 2 /fm 0. The model then explains the 
centrality dependence of S+U and Pb+Pb collisions at 
SPS energy. As the model parameters are fixed to repro- 
duce J/tp production in pA collisions, where deconfincd 
matter formation is unlikely, it was concluded that at 
SPS energy S+U/Pb+Pb collisions, J/tp's are absorbed 
in a nuclear medium 0, [l(| • 

If at RHIC energy, J/tp's are suppressed in a medium 
denser than the medium produced in SPS energy nuclear 
collisions, e 2 will increase. In a denser medium, the cc 
pair will interact more with the medium and per unit 
length will gain more square momentum. Parametric 
value of e 2 at RHIC energy then can indicate whether 
or not, a dense medium is produced. We note that 
due enhanced energy charmonium production increases 
at RHIC, but as shown in Figfl] that energy dependence 
is included in the model. 



III. J/tp SUPPRESSION IN CU+CU/AU+AU 
COLLISIONS AT RHIC 

PHENIX collaboration has measured the centrality de- 
pendence of J/tp suppression, in Cu+Cu and in Au+Au 
collisions, in two ranges of rapidity intervals, (i) —0.35 < 
y < 0.35 and (ii)1.2 < y < 2.2, |1| [27|. In FigM pre- 
liminary PHENIX data on the centrality dependence of 
nuclear modification factor (Raa) for J/tp, in Cu+Cu 
collisions are shown. Two rapidity ranges of data are 
not distinguished. We note that the present model is 
designed for central rapidity only. However, presently 
we ignore this limitation of the model. As it is evident 
from Fig|2l rapidity dependence of J/tp suppression is 
not large in Cu+Cu collisions. Using the CERN minimi- 
sation programme MINUIT, with e 2 as a parameter of 
the QCD based absorption model, we fit the data. The 
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FIG. 2: Preliminary PHENIX data on the centrality depen- 
dence of nuclear modification factor for J/tp in Cu+Cu col- 
lisions at RHIC. The solid lines is the best fit in the QCD 
based nuclear absorption model with e 2 =0.173 GeV 2 / fm. 



best fit is obtained with e 2 = 0.173 ± 0.007 GeV 2 / fm. 
The fit is shown in Figj2] (the solid line) . With the excep- 
tion of very peripheral collisions, PHENIX data on the 
centrality dependence of Raa, for J/tp in Cu+Cu colli- 
sions, are well explained the QCD based nuclear absorp- 
tion model. In very peripheral collisions, the model over 
predicts the suppression. Indeed, in very peripheral col- 
lisions, PHENIX data indicate enhancement rather than 
suppression of J/tp, presumably due to Cronin effect. 

As discussed earlier, if at RHIC energy, J/tp's are sup- 
pressed in a medium, denser than the medium created 
in SPS energy nuclear collisions, e 2 should increase. In 
Cu+Cu collisions, a contrary result is obtained. Com- 
pared to e 2 at SPS energy, at RHIC Cu+Cu collisions, 
e 2 decreases by a modest 7%. J/tp's are less suppressed. 
The result is consistent with the PHENIX measurement 
of J/tp suppression in d+Au collisions at RHIC [12]. At 
RHIC d+Au collisions, J/?/>-nucleon absorption cross- 
section, <Jj/ipN wl-3 mb, is less than the J/^-nucleon 
absorption cross-section at SPS energy, aj/^N ~4-5 mb 
[33l ]. Good fit to the centrality dependence of J/tp sup- 
pression in Cu+Cu collisions, data, in the QCD based 
nuclear absorption model, with e 2 close to the value at 
SPS energy, indicate that in Cu+Cu collisions, J/tp's are 
suppressed in a nuclear medium, much like the medium 
produced in nuclear collisions at SPS energy. 

Next we fit the PHENIX data on J/tp suppression in 
Au+Au collisions. Centrality dependence of the nuclear 
modification factor (Raa) for J/tp in the central rapidity 
region —0.35 < y < 0.35, is shown in Fig(3J Data points 
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FIG. 3: Preliminary PHENIX data on the centrality depen- 
dence of nuclear modification factor for J/ip in Au+Au colli- 
sions at RHIC. The solid line is the best fit in the QCD based 
nuclear absorption model with e 2 =0.146 GeV 2 / fm. 
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FIG. 4: Collision number dependence of the ratio of average 
number of gluon-nucleon collisions to the gluon-nucleon cross- 
section in Au+Au collisions. The solid and dashed lines are 
obtained with e 2 =0.14 and 0.18 GeV 2 / fm respectively. 



are few and error bars are also large. The best fit to the 
data is obtained with e 2 = 0.146 ± .014 GeV 2 / fm. It is 
shown as the solid line in Figj3l The fit to the data is 
not satisfactory. In very central collisions, the model over 
predict Raa and in mid central collisions the model under 
predict Raa- The best fitted value of e 2 = 0.146 ± .014 
GeV 2 1 fm is ~ 15% lower than the value required to ex- 
plain the PHENIX data on J/ip suppression in Cu+Cu 
collisions. Apparently, in Au+Au collisions, J/ip's are 
suppressed in a medium, less dense than that produced 
in Cu+Cu collisions. This is inconsistent with other re- 
sults at RHIC Au+Au collisions 0, E 01, S3]- As 
mentioned earlier, in RHIC Au+Au collisions, deconfincd 
matter can be formed. J/ip's can be suppressed in the 
deconfined matter. QCD based nuclear absorption model 
do not account for such a suppression. Unsatisfactory fit 
to the Au+Au data may be due to neglect of deconfincd 
medium production in Au+Au collisions. 



IV. p T BROADENING OF J/ip IN CU+CU AND 
AU+AU COLLISIONS 

It is well known that in pA and AA collisions, the sec- 
ondary hadrons generally shows a px broadening [34], [35| . 
Pt broadening of J/i/j's, in S+U and Pb+Pb collisions 
are well explained in the QCD based nuclear absorption 
model. Recently PHENIX collaboration has measured 
Pt broadening of J/ip in Cu+Cu and Au+Au collisions 
[2& It is interesting to compare the QCD based 



nuclear absorption model predictions with the PHENIX 
data. 

The natural basis for the pr broadening is the initial 
state parton scatterings. For J/ip's, gluon fusion being 
the dominant mechanism for cc production, initial state 
scattering of the projectile/target gluons with the tar- 
get/projectile nucleons causes the intrinsic momentum 
broadening of the gluons, which is reflected in the pr 
distribution of the resulting J/ip's. Parameterising the 
intrinsic transverse momentum of a gluon, inside a nu- 
cleon as, 

f{q T ) ~ exp(-q T / < q 2 T >) (9) 

momentum distribution of the resulting J/ip in NN col- 
lision is obtained by convoluting two such distributions, 

f J N f N (Pt) ~ exp(-p 2 T / < p\ >iit) (10) 

where < p\ > J Jj%=< q T > + < q\ >■ In nucleus- 
nucleus collisions at impact parameter b, if before fu- 
sion, a gluon undergo random walk and suffer num- 
ber of subcollisions, its square momentum will increase 
to q 2 ^ — * q T + A^o, <5o being the average broadening in 
each subcollisions. Square momentum of J/ip then easily 
obtained as, 

< Pt At (b) =< P% > J Jt +SoN AB (h) (11) 



5 



where Nab (b) is the number of subcollisions suffered by 
the projectile and target gluons with the target and pro- 
jectile nucleons respectively. 

Average number of collisions Nab (b) can be obtained 
in a Glauber model [351 ] - At impact parameter b, the po- 
sitions (s, z) and (b — s, z') specifies the formation point 
of cc in the two nuclei, with s in the transverse plane and 
z, z' along the beam axis. The number of collisions, prior 
to cc pair formation, can be written as, 



fit is obtained with a gN S Q = 0.03 ± 0.51 GeV 2 . RHIC 
data do not show any evidence of pt broadening, as in- 
dicated by very small value of cr g jv<5o- In Fig(5j model 
predictions with the central value, <7 s at<5o =0.03 GeV 2 
are shown. The solid and dashed lines are for Au+Au 
and Cu+Cu collisions respectively. The predictions for 
Cu+Cu collisions closely match that for Au+Au colli- 
sions and cannot be distinguished. 



N(b,s,z,z')= <JgN dz A pA(s,z A ) (12) 

J — oo 

+<r gN / dz B pB(b - s,z') 

where a g N is the gluon-nucleon cross-section. Above ex- 
pression should be averaged over all positions of cc for- 
mation with a weight given by the product of nuclear 
densities and survival probabilities 5*, 

/p oo />oo 
d 2 s dzp A {s,z) dz'p B (b- s,z') x 
J -co J —oo 

S{b,s)N(b,s,z,z')/ I ' d 2 s f dzp A (s,z)x 

J J — oc 

dz' p B (b — s, z')S(b, s, z, z') (13) 



In FiglU the centrality dependence of the ratio 
Nab I VgN i in Au+Au collisions are shown. The solid 
and dashed lines corresponds to e 2 =0.14 GeV/fm and 
0.18 GeV/fm respectively. NAB/<^gN do not show large 
dependence on e 2 . Even though e 2 differ by 25%, Nab 
differ by less than 3% in central collisions. In less central 
collisions, the difference is even less, pt broadening of 
J '/if> will not depend much on the exact value of e 2 . 

Pt broadening of J/f/Vs in AA collisions depends on 

two parameters, (i) < p\ >jvjv> the mean squared trans- 
verse momentum in NN collisions and (ii) the prod- 
uct of the gluon-nucleon cross-section and the average 
parton momentum broadening per collision, a g N^o- < 
Pt > 'nn i s a measured in RHIC energy p+p collisions, 



< p 2 T > J J%= 4.2+0.7 GeV 2 . The other parameter, a. 



gN°0 



is essentially non-measurable, as gluons are not free. Its 
value can only be obtained from experimental data on 
Pt broadening of J/ij). At SPS energy S+U/Pb+Pb col- 
lisions CT 9 Ar^o is estimated as 0.442 ± 0.056 GeV 2 [13]. 
PHENIX data on J/tjj pt broadening can be used to es- 
timate its value at RHIC energy. 

In Fig[5l PHENIX data on the centrality dependence 
of mean square transverse momentum < p T >, in Cu+Cu 
and in Au+Au collisions are shown. For comparison pur- 
pose, < p T > in in p+p and in d+Au collisions are also 
shown. As data points are few, we do not fit the individ- 
ual Cu+Cu or Au+Au data sets. Rather we fit the com- 
bined data sets. We fix < p T >nn at the measured cen- 
tral value, < p\ >nn — 4.2 GeV 2 , and vary ct 9 jv^o- Best 
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FIG. 5: J/ip mean square transverse momentum as a func- 
tion of collision number, in p+p, d+Au, Cu+Cu and Au+Au 
collisions are shown. The solid and dashed lines are fit to the 
Au+Au and Cu+Cu data respectively. 



V. SUMMARY AND CONCLUSIONS 

To summarize, we have analysed the (preliminary) 
PHENIX data [24 l25| on the centrality dependence of 
J/ip suppression and pt broadening, in Cu+Cu and in 
Au+Au collisions at RHIC energy, y/s=200 GeV. The 
data are analysed in the QCD based nuclear absorption 
model d, [13, LU], [l3,G3. I n t ne m odel, J/^'s suppression 
is controlled by a parameter e 2 . Larger the e 2 , more is 
the suppression. Centrality dependence of Jji\) suppres- 
sion at SPS energy require £| PS =0.187 GeV 2 j jm. With 
e 2 as a parameter, we have fitted the PHENIX data on 
the centrality dependence of J/ip suppression in Cu+Cu 
and in Au+Au collisions. Cu+Cu data are well explained 
with e 2 = 0.173 ± 0.007 GeV 2 /fm, close to the SPS en- 
ergy value. In Cu+Cu collisions, J/ip's are suppressed in 
a medium much like the medium created in SPS energy 
collisions. No exotic, high-density matter is created in 
Cu+Cu collisions. 

Centrality dependence of J ftp suppression, in Au+Au 
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collisions, is not well explained in the model. Best fitted 
value, e 2 = 0.146 ± .014 GeV 2 /fm, over predict the sup- 
pression in mid-central collisions and under predict the 
suppression in very central collisions. We conclude that 
in Au+Au collisions, J/ip's are suppressed in medium un- 



like the medium created in SPS energy nuclear collisions 
or in Cu+Cu collisions at RHIC energy. We have also 
analysed the PHENIX data on pt broadening of J/ ^ in 
Cu+Cu and Au+Au collisions. RHIC data onpx broad- 
ening do not show any evidence of pr broadening. 
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